Field Theory of the Random Flux Model 
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The long-range properties of the random flux model (lattice fermions hopping under the influence 
of maximally random link disorder) are shown to be described by a supersymmetric field theory of 
non-linear o model type, where the group GL(n|n) is the global invariant manifold. An extension 
to non-abelian generalizations of this model identifies connections to lattice QCD, Dirac fermions 
in a random gauge potential, and stochastic non-Hermitian operators. 



Quantum disordered systems are typically realised in 
Hamiltonians of the general form H = Hq + V, where Ho 
models the underlying "clean" system, and disorder is in- 
troduced via the randomly distributed Hermitian opera- 
tor V. Sometimes, however, it is preferable to implement 
disorder in terms of unitary stochastic operators and to 
consider Hamiltonians of the type 



1 U%j c j > 



(i) 



where (ij) denote neighbouring sites of a d-dimensional 
hypercubic lattice, the c's represent A-component lattice 
fermions, and Uij represent A-dimensional unitary ma- 
trices residing on the links of the lattice. Stochasticity 
is introduced by drawing the ?7's from a random dis- 
tribution (albeit subject to the Hermiticity requirement 



U, 
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Hamiltonians of the type (|l|) are commonly 
referred to as random flux (RF) Models, a denotation we 
will also adopt for the cases JV^ 1. 

RF-models appear in a variety of different contexts: 
The 2d A = 1 version describes the dynamics of lattice 
fermions subject to a random magnetic field or, more ac- 
curately, a random vector potential Jl] ||] . As well as the 
gauge theory of high T c superconducitvity , this model 
has been discussed in connection with the physics of the 
half- filled fractional quantum Hall phase , as well as the 
spin-split Landau level |3|,^| . Identifying the two fermion 
components of the A = 2 RF-model with a spin degree 
of freedom, (|l|) describes the propagation of lattice elec- 
trons on a spin-disordered background, a situation that 



occurs, e.g., in connection with the physics of manganese 
oxides Identifying the three fermion components of 
the A = 3 model with a color degree of freedom, 



represents a prototype |10| of the strong coupling lattice 
QCD-Hamiltonian. 

Superficially, (|l|) appears to fall into the general class 
of (bond disordered) Anderson Hamiltonians. That con- 
jecture indeed holds true provided one stays away from 
the middle of the tight-binding band, e = 0. Upon ap- 
proaching e = 0, however, the phenomenology of the RF- 
models begins to differ drastically from a conventional 
disordered fermion systems. In spite of intensive numeri- 
cal and analytical investigation , central aspects of 
these deviations are not yet fully understood. For exam- 
ple, the key question of whether or not the 2d RF-model 
possesses a band center extended metallic phase has not 
yet been settled; apart from the fact that the average 
density of states (DoS) diverges upon approaching e = 0, 
much of the structure of even that basic observable re- 
mains unknown. 

The purpose of this Letter is 2-fold: Firstly we provide 
new information regarding the band center behavior of 
the RF-model. Secondly we wish to discuss a diverse 
network of interconnections that exist between the RF- 
problem and related areas of current research interest. 

Both aspects of that program are based on the result 
that the long-range behavior of average n-point Green 
functions, (G ± (e + e 1 )...G ± {e + e„)), of the RF-model 
can be obtained from a supersymmetric field theory de- 
fined by the action 



S[T] 



d str(0r _1 0T) + ic 2 stT (e(T + T" 1 )) + c 3 ( str (T^dT))' 



■S h [T] 



(2) 



where T € GL(n\n) (the group of invertible supermatri- 
ces of dimension 2n), 'str' is the standard supertrace, and 
the matrix e = diag(ei, . . . , e n ). The contribution Sb JTl] ] 
represents a boundary action that depends on the values 
of the fields T at the corner points of the lattice. 

Eq. (^) is derived under the assumption of maximal 
unitary randomness, i.e. all Uij £ U(A), independently 



distributed according to the Haar measure. In this case, 
the constants a = Na 2 ~ d /M, c 2 = N(2d - \) 1 l 2 aT d j<±d, 
C3 = a 2 ~ d C /16d where a represents the lattice spacing, 
and C denotes a geometry-dependent numerical constant 
0(1). Below we will argue that the structure of the field 
theory is actually disorder independent, i.e. that RF- 
models are generally described by (§) p2[ , where the 
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strength of the disorder manifests itself merely in the 
value of the coupling constants. 

Below we will outline how, starting from the 'micro- 
scopic' Hamiltonian ([!]), the effective description (^) is 
derived. However, before turning to that more technical 
part of the discussion, we first address the question of 
what kind of information can be gained from the field 
theory. Our main goal will be to demonstrate that the 
action (Q) represents a quantitative implementation of 
the network of connections displayed in Fig. 1. By ex- 
ploring different links, we will discuss some characteristic 
features of the field theory. 



Non-Hermitian QM 



Random Dirac Fermions 
(Gauge Disorder) 



Lattice 
QCD 



RF-Model 



Sublattice 
Models 



Chiral RMT 



FIG. 1. Connection between the RF-model and various re- 
lated systems governed by the presence of chiral symmetries. 
The low energy limit of all models is universally described by 
Chiral Random Matrix theory. 

Chiral Random Matrix Theory (ChRMT): As usual 
with field theories of disordered systems, the low energy 
regime of (g) (energies e < ci/(c 2 L 2 )) is governed by 
spatially constant field configurations T = const., 

So[T ] = -^str (e(T + T^ 1 )) + S h [T ], (3) 

where po is the bulk mean DoS of the system. Corre- 
lation functions computed with respect to the first con- 
tribution to the action (^) coincide with those otherwise 
obtained for the chiral unitary random matrix ensemble 
ChGUE , (the ensemble of symmetry AIII in the 

classification scheme of Ref. |17|]), i.e. the ensemble of 
block off-diagonal matrices 



At 



(4) 



where A is complex random. In particular, the mean 
DoS is found to vanish as e — * on a scale set by the 
mean level spacing. The connection to ChRMT follows 
readily from the fact that the RF-Hamiltonian possesses 
a chiral symmetry: Partitioned into two nested sublat- 
tices A and B, the bipartite lattice Hamiltonian assumes 
a block off-diagonal form (^) in an A/£?-decomposition. 
To the best of our knowledge, the ramifications of the 
chiral structure on the physical properties of RF-models 



was first reported in Ref. M. In passing we note that, 
in contrast to conventional disordered systems, the low 
energy limit of the RF-model is not absolutely universal: 
The fine structure of the DoS close to e — depends on 
the 'parity' of the lattice, i.e. on whether the number of 
sites is even or odd p8[ . Without going into details, we 
remark that the information about this effect is encoded 
in the boundary term S\> . 

Non-Hermitian Operators: To investigate problems in- 
volving non-Hcrmitian stochastic operators A ^ A' one 
commonly introduces an operator like the one shown 
in (0), i.e. an Hermitian auxiliary operator of twice 
the dimension of the original problem Put dif- 

ferently, non-Hermitian Hamiltonians possess an inbuilt 
chiral structure implying that their low-energy univer- 
sal properties coincide with those of manifestly chiral 
problems like the RF-model. That in turn means that 
the basic structure of the low energy field theory @) of 
the RF-model (a system with broken time-reversal in- 
variance) should coincide with that of the (time-reversal 
non-invariant version of the) non-linear cr-model of non- 
Hermitian problems introduced in Ref. pjj |. To make 
that connection explicit, we introduce the auxiliary ma- 
trix variable Q = exp(VFeri/2)(s Cs5 03) exp(— Wa\/2), 
where W = InT, Oi are Pauli matrices operating in the 
block space of (|j), and s — diag(sgnlmei, . . . , sgnlme„). 
One may check by direct comparison that for n = 1 (the 
case there considered), the matrices Q are equivalent to 
the degrees of freedom employed in Ref. |2l[] . When rep- 
resented in terms of Q's, (^|) assumes a form similar to 
a standard [f22[ non-linear cr-model, albeit one of novel 
symmetry [eTT That the connection is not incidental, 
but rather extends to the more complex variants n > 1, 
follows from a) the above mentioned fact that both non- 
Hermitian and RF-type problems possess a chiral struc- 
ture, and b) that only three fundamentally different 0- 
models with chiral symmetry (corresponding to the cases 
of broken time-reversal invariance, broken spin rotation 
invariance, and invariance under both operations) exist. 
For a more thorough discussion of these symmetry as- 
pects, we refer the reader to the original Ref. O] . 

Weakly Disordered Sublattice Models: Leaving the ran- 
dom matrix regime and turning to the more complicated 
spatially extended problem, it is important to notice that 
the field theory (||) has a closely related precursor: Analy- 
sis of a weakly disordered sublattice model led Gade |23j | 
to a boson replica version of the present model, i.e. a 
theory over fields T E GL(nR)/\J(nR), where R — > is 
the number of replicas. The action for these fields coin- 
cides with (|^) , save for the absence of the boundary term, 
and the important difference that, due to the weakness of 
the disorder, the coupling constant c\ was parametrically 
larger than one. 

Various conclusions concerning the physical behavior 
of the RF-Hamiltonian, most notably about its localiza- 
tion behavior, can be inferred directly from Ref. [E3l. 
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There it was shown that the conductance of the weakly 
disordered 2d model at the band center (which is es- 
sentially determined by the coupling constant c%) did 
not change under one-loop perturbative renormaliza- 
tion. This observation suggests that a non-localized state 
might exist in the middle of the band. Since the stabil- 
ity of the perturbative RG relies merely on the small- 
ness of the parameters l/c\,c^/c\ -C 1, its results can be 
straightforwardly carried over to the N ^ 1 non-abelian 
RF model: The one-loop renormalization indicates that, 
at least for N ^> 1, the strongly disordered RF model 
exhibits metallic behavior at the band center. (It is in- 
teresting to note that, according to the connections sum- 
marized above, the unusual localization properties of the 
zero energy states of the RF-model characterize those of 
all eigenstates of a stochastic non-Hermitian operator.) 

It was also predicted in Ref. (23) that the 2d DoS di- 
verges upon approaching the middle of the band. In or- 
der to understand on which energy scale that divergence 
sets in, and how it will eventually be cut off deep within 
the random matrix regime, one would have to superim- 
pose perturbative RG techniques onto a non-perturbative 
treatment of the low-energy regime, a task that is well 
beyond the scope of the present paper. 

Finally we notice that, in an RG-sense, finite energies 
Ci represent a relevant perturbation. Renormalization 
of the field theory leads to a flow away from the chiral 
band center limit eventually leading to the standard uni- 
tary universality class. This result is consistent with the 
analysis of Ref. where it was shown that continuum 
fcrmions (i.e. the analogue of lattice fermions close to 
the bottom of the band) subject to a weak random field 
map onto a unitary cr-model. 

The tendency of sublattice models to exhibit band cen- 
ter delocalized behavior persists even in the (quasi) Id 
case: It was shown in Refs. ||^4| that Id sublattice mod- 
els with N even exhibit conventional localization behav- 
ior whilst for N odd a delocalized mode remains in the 
band center. This parity effect is closely related to the 
odd/even phenomenon mentioned above in connection 
with the mean DoS, and indeed it is the boundary ac- 
tion that is responsible for the quasi Id derealization 
phenomenon within the cr-model formulation. 

Lattice QCD and Random Dirac Fermions: Besides 
Gade's model, the field theory (||) has at least two other 
close relatives: In QCD, (||) has been suggested on phe- 
nomenological grounds as relevant for the determination 
of the low energy spectrum of the Dirac operator [ B5[ . 
In that context, the base manifold is 4 + 1-dimensional 
whilst the fields T € U(nf + 1|1), where n/ is the num- 
ber of quark flavors. For a comprehensive discussion of 
the QCD-analogue of (|J), its connection with ChRMT 
and its relevance for lattice QCD-analyses, we refer the 
reader to Ref. p6|. The similarity between the theories 



is again a manifestation of the universality of chiral a- 
models or, more physically, the universal consequences 
chiral symmetries have for the long-range properties of 
random systems. (In QCD, 'randomness' is represented 
by gauge field fluctuations in the Yang-Mills Hamilto- 
nian.) 

Secondly, an analogue of the action (0) (again without 
the boundary operator oj) with a Wess-Zumino-Novikov- 
Witten WZWN term 

-c 5 f dCe^str (Td( i T~ 1 Td tJ T~ 1 Td v T~ 1 ) 
Jo 

was obtained from a lattice model of random Dirac 
fermions ^7j. The connection to the RF-model follows 
from the fact that the clean limit of the Dirac model can 
effectively be described in terms of a model of 7r-flux lat- 
tice fermions. (The Dirac-structure of the clean 7r-flux 
model is also responsible for the occurrence of a WZMN 
operator.) 

Having reviewed some essential features of the field 
theory (^|), we finally outline how it is obtained from 
(|I|). That in this Letter, the derivation is not formu- 
lated in more detail is motivated by the observation that 
not only the degrees of freedom but also the structure 
of the field theory is, to a large extent, dictated by as- 
pects of symmetry: By analogy to the situation for the 
'conventional' supersymmetric c-models | j22f , there are 
only a few GL(n|n) invariant operators with < 2 gradi- 
ents (namely the ones appearing in ^ plus the WZWN 
operator |28|). Thus, the 'only' job that is left for a mi- 
croscopic derivation is to decide whether the operators 
permitted by symmetry are actually realized in the field 
theory, and to fix their coupling constants. Here we re- 
strict ourselves to a brief outline of that analysis. Details 
of the calculation will be presented in a separate publi- 
cation. 

(i) As usual in the construction of field theories of dis- 
ordered problems, we first represent Green functions of 
the problem in terms of a Gaussian integral over a field ?p. 
Choosing supersymmetry as a way to normalize the re- 
sulting functional integrals to unity, the first step exactly 
parallels the constructions reviewed in Ref. (^|| . (ii) Next 
we average over the set {Uij}. At that stage significant 
deviations from the standard treatment of Hermitian dis- 
ordered operators occur. A method of exactly averaging 
over (extended) models involving unitary stochasticity 
has been introduced in Ref. |2{],|3(| and christened the 
'color-flavor transformation'. Following that reference 
we eliminate the disorder at the expense of introduc- 
ing a pair of auxiliary fields {(Zij, Zij)} (which play a 
role analogous to the Hubbard- Stratonovich field Q com- 
monly employed in a- model constructions), (iii) Inte- 
grating out the ^'s we are left with the action 
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s[z,z} = ~nj: 



(i£A,j£B) 



str In (1 



■N J2 str m U ■ 



E 



Zi 



+ N Yl str In (e 

jeB v 



(5) 



where the notation j G Nj indicates that j is summed 
over all nearest neighbours of i. (iv) Subjecting (||) to 
a saddle-point analysis, the fields are conveniently pa- 
rameterized as (Z,Z) = (ixPT, ia;T _1 P), where a; is a 
constant, and T,P G GL(n|n) respectively have the sig- 
nificance of Goldstone, massive modes of the theory, (v) 
Integrating out P we find that, in contrast to standard 
cr-model analyses (and in accord with the construction of 
Gade's action |25(] ), a residual coupling between massive 
and Goldstone modes exists; it gives rise to the C3-term 
in (g). (vi) The remaining, pure Goldstone action is sub- 
jected to a gradient expansion which results in (|2|). 

Summarizing, we have derived an effective field theory 
for the maximally disordered RF-model. The theory has 
a status analogous to the supersymmetric non-linear a- 
models of 'conventional' disordered Fermi systems, but 
its behavior is substantially different, a fact that is read- 
ily traced back to the presence of a chiral symmetry. It 
was shown that the formalism provides a platform from 
which interconnections to a variety of other recently in- 
vestigated chiral problems can be conveniently analysed. 
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